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Magnonics is an emerging field with potential applications in classical and quantum information 
processing. Freely propagating magnons in two-dimensional media are subject to dispersion, which 
limits their effective range and utility as information carriers. We show the design of a confining 
magnonic waveguide created by two surface current carrying wires placed above a spin-sheet, which 
can be used as a primitive for reconfigurable magnonic circuitry. We theoretically demonstrate 
the ability of such guides to counter the transverse dispersion of the magnon in a spin-sheet, thus 
extending the range of the magnon. A design of a magnonic directional coupler and controllable 
Michelson interferometer is shown, demonstrating its utility for information processing tasks. 

PACS numbers: 75.10.Pq, 03.67.Hk, 05.60.Gg, 75.30.Ds 


I. INTRODUCTION 

During the past two decades, magnonics has be¬ 
come a major field of theoretical and experimental 
researchP^. The potential use of spin waves in several 
technologies, like computing (quantunPand classicaP^^, 
communicatiorP and caloritronicP has been investigated 
thoroughly. Magnon-magnon scattering also provides an 
interesting test bed f or exp loring particle-particle inter¬ 
actions more generall} hO | ii |, The core idea of magnonics 
is to manipulate spin waves in ferromagnetic or antifer¬ 
romagnetic materials, and i n doing so to produce devices 
with different functionalitieP^*^. This manipulation is 
usually achieved by the application of external fieldP^®^, 
by crafting the material itself (remo ving th e material on 
sides the achiev e a p hysical channel (pHlSl or some com¬ 
bination of botffi^Ilffl 

In a magnonic device there is no net transport of charge 
as opposed to electronic devices, where current flow leads 
to inevitable heating. For this reason, magnonic technol¬ 
ogy promises the advantage of reduced power consump¬ 
tion compared to its electronic counterparts.^^ 

A spin wave is a propagating magnetic distur bance 
in an otherwise perfectly ordered magnetic mediunP^^. 
When quantised, the associated particle is a magnon. 
Here we are concerned with the propagation of magnons 
confined to a two dimensional array of spin-1/2 particles 
(spin sheet). Within the sheet, each spin interacts with 
its nearest neighbour through exchange or dipole-dipole 
coupling. Due to this interaction they either align (ferro¬ 
magnetic coupling) or anti-align (anti-ferromagnetic cou¬ 
pling) themselves to their nearest neighbours, resulting 
in a global order in spin direction as the lowest energy 
state. Any local disturbance in the global order creates 
propagating magnons in the sheeif^^. 

Freely propagating magnons are subject to damping 
and dispersion, which limits the distances over which the 
magnon can propagatd^^HSSI The effects of dispersion 
can be mitigated by confining the magnon in the direc¬ 
tion transverse to the travel direction. Typically, such 


confinement is achieved by lithographically removing ma¬ 
terial to leave a magnonic channeF^*^, which minimises 
the transversal dispersion. However, such approaches are 
not suited to reconfigurable devices due to the permanent 
structural changes induced in the material. 


(a) 


PMC. 





(b) 



FIG. 1. (a) Schematic of spin guide. Black lines represent cur¬ 
rent carrying wires with current flowing in opposite direction. 
Black translucent surfaces are magnitude field contours. The 
surface wires define a confining potential in the direction of 
the wire, (b) Schematic of the magnonic Michelson interfer¬ 
ometer. The black layer shows a two-dimensional ferromag¬ 
netic sheet on a substrate. Brown lines are current carrying 
wires and the arrows show the direction of the current in the 
wires. One arm has a dynamic magnonic crystal (DMC) at 
its end for the purpose of phase control. 
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Here we describe a method of confining and transport¬ 
ing a magnon in a two-dimensional ferromagnetic film, 
which does not require any structural modification of the 
film itself. In our model, the confinement is achieved by 
the temporally invariant magnetic field from two paral¬ 
lel current-carrying wires placed on the surface at some 
distance from each other. Schematics of such a device is 
shown in fig. (a), where fig. (b) shows the design of 
magnonic Michelson Interferometer. 

When the current direction in the wires is opposite, the 
resulting magnetic field forms a potential well, creating 
a channel for the magnons, which is termed a spin-guide. 
Previously we have shown analogous confinement in one¬ 
dimensional spin chains using a surface-gate array, where 
the transport was achie ved th rough adiabatic temporal 
variation of the potentiaP^^^. 

Our results enable reconfigurable magnonic circuitry 
of several, potentially interacting, magnon channels that 
preserve phase and thereby allow quantum information 
processing. The channels can be reconfigured by adapt¬ 
ing the current in the wires that create the channel. One 
can consider multiple spin-guides on a single surface that 
can be switched on or off as desired. This ability can be 
used to build complex magnonic circuitry for all-magnon 
information processing. More generally, our work is mo¬ 
tivated by efforts to demonstrate traditional quantum op¬ 
tical effects in non-traditional tight-binding system^^llESl 

To illustrate the utility of our scheme, we present a full 
design for a guided magnon version of a Michelson inter¬ 
ferometer. The Michelson interferometer is an important 
element for high-precision sensing and can also be used 
as a switching primitive for all-magnon logic gate^^. 

This paper is set out as follows. We first introduce the 
Heisenberg Hamiltonian in two dimensions and a magnon 
confining potential. We then consider the extension to 
a directional coupler and finally a magnonic Michelson 
Interferometer. 


II. FREE AND GUIDED MAGNON 
PROPAGATION IN A 2D HEISENBERG SHEET 

To calculate the magnon propagation, we introduce the 
Heisenberg Hamiltonian with nearest neighbour interac¬ 
tion in a two dimensional square lattice, with h = \. 

H = —J ^ ^ [Sij-.Si+ij- -|- ^ ^ (1) 

i,j i,j 

S = ,a-^) is the total spin operator, where 

cr“,errand are the Pauli spin operators. J is the 

strength of nearest neighbour interaction. The nature 
of the interaction can be exchange or dipole-dipole. If 
both types of couplings are present then J becomes the 
effective coupling with the magnitude equal to the sum 
of both coupling, i.e. J = Jex + Jd, Jex is the exchange 
coupling and Jd is the dipole-dipole couplin^^. 

The third term in the Hamiltonian [eq. 0 ] is the on¬ 
site energy of (i, j)th spin due to the local external mag¬ 


netic field , with where 7 is the gy- 

romagnetic ratio. 

We do not take into account the second-nearest neigh¬ 
bours i.e. (f, j) to {i -I- l,j -I- 1) coupling. For magnons 
propagating close to parallel to either of the principle 
axes, the effects of such couplings can be approximated 
by an effective increase of the nearest-neighbour cou¬ 
pling. Through an analytical derivation (see Appendix) 
we demonstrate that the effect of these interactions is 
equivalent to that of having a stronger J coupling, for 
magnons travelling in either x or j/-direction. For such a 
case, to account for second-nearest-neighbour coupling, 
we substitute J = Je + ‘2Jd in eq. Q, where Je is the 
standard edge coupling and Jd is the diagonal coupling 
on a 2D square lattice. 

The dispersion relation of a freely propagating magnon 
in a two dimensional sheet is 

^k^,ky = = 4J - 2J[cos{kxa) -f cos(fcya)], (2) 

where ujk^,ky is the angular frequency of a particular 
eigenmode, kx and ky are the components of iT, wave 
numbers along x and y directions and a is the spin¬ 
spacing. The dispersion curve of a square lattice along 
the high symmetry lines is shown in fig. [|:a,b). 

To illustrate the free magnon dispersion we considered 
a magnon with the initial wave function 


i’ixi,yj) =exp 


(a:o - XiY 

^^x 


{Vo-Vjf 


eyi]i{-ikxXi - ikyyj). 


( 3 ) 


Equation (|^ corresponds to a two-dimensional Gaussian 
wavepacket propagating in the x and y-direction with 
the initial group velocity corresponding to wave vector 
K = kx X + kyjj. (j)x and (py are the spatial standard 
deviations in the x and y-directions and (xq, yo) is the 
center of the wave packet. 

In an experimental setup, the standard way to excite 
magnons in a thin film is by the microwave induction 
techniqu ci^ l ^^i i^. A surface antenna is subjected to AC 
current and the Oersterd field of the antenna introduces 
propagating spinwaves inside the sheet. The frequency 
and the wavelength of the excitation can be controlled di¬ 
rectly by the frequency of the AC current. When such an 
antenna is used on the top of the guide to excite magnons, 
only the confined modes will travel along the guide and 
all unconfined modes will leak out and disperse away. 
This approach can be thought of as being analogous to 
coupling into a few-mode optical fiber. 

The required iT-vector to achieve a certain group ve¬ 
locity of the magnon can be calculated by the following 
relation 


Vg =VgX -f Vgi) = dk,UJk,,kyX + ,/c„ J/ 

=2J a[si'a(kxa)x + sin(fcya)y]. 


( 4 ) 
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(b) 





FIG. 2. (a) Brillouin zone of a square lattice with the lines of high symmetry marked in bold, (b) Dispersion relation of the 
square lattice along high symmetry lines. Red and black ellipses show the areas of high and low dispersion, respectively, (c) 
Overlaid snap shots of a propagating magnon in a square-lattice, while the color axis is I'i/'p. Rapid transversal dispersion is 
apparent where the longitudinal dispersion is not noticeable at this distance. This is because the magnon was chosen with a y- 
momentum, which is in a low dispersion region, as marked in black in (b). (d) Magnon confined by the surface current-carrying 
wires for the same initial state as (c). Red arrows show the direction of the current in the wires and white lines represent wires 
(d = 20a, Wg = 20a and emin/J = 0.1). After the confinement there is no transverse spreading of the magnon. 


In our simulation we initiated the magnon with a y- 
velocity Vg = 2aJy, which corresponds to kyU = 7r/2 
and kx = 0. 

Figure [^c) shows the evolution of a free magnon along 
the spin sheet. The time evolution was achieved through 
solving the discrete version of the Schrodinger equation 
and hard wall boundary conditions were implemented. 
We observe strong dispersion in the transverse direction 
and weak dispersion in the lateral direction. This is due 
to the fact that the two directions correspond to the dif¬ 
ferent parts of the dispersion curve. In the propagation 
direction, the group velocity is closer to the linear part 
[fig- l^b) black ellipse] of the dispersion curve, hence it 
shows reduced dispersion, whereas the zero momentum 
region is a high dispersion region [fig. Hb) red ellipse] 
and hence spreading is more rapid. 

For our confinement scheme, we use two current carry¬ 


ing wires with equal but antiparallel currents, placed at 
some distance above the plane of the spin-sheet, as shown 
in fig. i a). Figure!^ (d) shows the evolution of a magnon 
inside the guide. As expected, there is negligible trans¬ 
verse dispersion, however the longitudinal dispersion is 
unchanged. The white lines show the wires and the red 
arrows show the direction of current in the wires. The 
longitudinal initial state for the confined case were same 
as of the free case. The transversal profile was chosen to 
be the ground state of the confining potential. 

In addition to the field of the wires, we also considered 
a global z-field, which is large compared to the wires’ 
magnetic field. This large field energetically separates the 
total ground state from the single-excitation subspace, 
and also allows the secular approximation. The secular 
approximation allows us to neglect all terms from non-z 
components of the wire-field because they are averaged 
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FIG. 3. (a) Cross section of a spin-guide. Two circles show the position of the wires, where the inner circle and cross inside 
each circle signifies the direction of the current. The horizontal line is the plane of the spin sheet. 2wg is the separation between 
wires and d is the distance between the spin sheet and the wires, (b) Onsite energy of spins inside the sheet due to individual 
( — ■— and — •— ) and combined ( — C— ) magnetic field of two current carrying wires. Two vertical (dotted) lines show 
the position of the wires. The two wires form a potential well with depth £min, which can guide magnons.(c) A pseudocolor 
plot of the depth, emin, due to the magnetic field of the wires, located at the center of the potential. The black region is the 
region where the potential starts to split into two separate potentials and therefore leads to more complicated dynamics. The 
inset shows the shape of such a split potential with the parameters that lie inside the black region(ci = 80a, tUg = 400a). The 
white line marks the position of the maximum potential depth. 


out by the strong z-field. Hence we perform all calcula¬ 
tions in the single excitation subspace and ignore the x 
and y-magnetic field components due to the wires. The 
functional form of the combined z-field of the wires at 
some point x inside the sheet is 


Bz{x) 


/xq/ x-Wg 

27r i?2 


X -\-w„ 
R'2 


( 5 ) 


where R = ^/d'^ + {x — WgY, R' = (P + {x -\- WgY, d 
is the distance between the sheet plane and wires, Wg is 
the half-width of the guide (half-separation between the 
wires), /io is the permeability of free space and I is the 
current in the wires. Figure [^a) shows a cross section 
of the sheet and wires. Two circles represent the wires 
running perpendicular to the plane of the page and the 
horizontal line is the plane of the spin-sheet. The current 
direction in each wire is opposite to each other. 

Figure [^b) shows the z-component of the magnetic 
field inside the sheet due to the individual wires and their 
combination. The maximum depth occurs at the middle 
of the wires. Figurej^c) shows the depth of the potential 
(Smin) as a function of Wg and d. The white line marks 
the maximum depth of the potential and hence defines 


an optimal choice of Wg and d, which is Wg = d, i.e. the 
distance of the wires from the spin sheet should be the 
half-width of the guide to achieve the maximum poten¬ 
tial depth, Emin- If the half-width, Wg, of the guide is 
significantly larger than its distance from the spin sheet, 
d, then the potential splits into two separate potentials 
(inset fig. [^. The region of split potential is marked in 
black in fig. ib). On the other hand if the spin-sheet to 
guide distance, d, is much larger then its half-width, Wg, 
this results in a shallow potential from which the magnon 
can escape more easily. Therefore, we restrict ourself to 
the choice of geometry where Wg = d. Substituting this 
and X = 0 in eq. ([^, we get an expression for the poten¬ 
tial depth = /x//27rd. The energy of a spin due to 

the which corresponds to the centre of the guide, 

see fig. |^(a), is Emin = = "fh^I/{2TTd) = I/{Kd), 

where 7 is the gyromagnetic ratio, /x is the magnetic per¬ 
meability of the material and k = 2tt j ('jHfi). This expres¬ 
sion gives a clear connection between geometry, energy 
and the current in the wires. 

For the rest of the paper we will use the potential depth 
£min as the independent variable. For any given scale 
d = Wg, the potential depth is proportional to current 
I. The required Emin to achieve a certain confinement 
profile scales linearly with the strength of the J-coupling 
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FIG. 4. (a) Confinement factor as a function of £min[J] for several values of d. To achieve the same CF, a potential with smaller 
d requires larger £min[J]. (b) Number of confined modes inside the potential. The criterion of confinement was chosen to be 
the value of CF. When CF of a particular mode reaches 0.9, we considered it a confined mode. A similar dependence on d 
holds true here, where smaller d needs a stronger potential to confine the same number of modes as a guide with larger d. (c) 
Ten lowest eigenvalues as a function of £min[J], d = Wg = 150a. A value below zero is a bound mode. As the potential depth 
increases, more and more modes become bound modes, (d) The difference between ground state and the first excited state. 
The difference also increases with the increasing depth of the potential. 


and therefore Emin is represented in units of J. With the 
knowledge of J of a particular system and the geometry 
of the guide, one can calculate the current that is required 
in the wires to produce the required Emin- 

Our spin guide model provides confinement analogous 
to that in optical wave guides and it can show similar 
functionalities, for example, bending and coupling. Al¬ 
though the shape will have an effect on the properties 
like adiabatic bending and coupling, which we discuss in 
detail in the next section, the exact shape of the confin¬ 
ing potential is not important for the guiding. As long 
as the potential is strong enough to create one or more 
confined modes, guiding can be achieved. 


III. CONFINEMENT FACTOR AND BEND 
LOSS 

Confinement factor is a commonly used term in semi¬ 
conductor laser physics, which is defined as the ratio of 
the electric field in the active region to the total electric 
fielcP. We use the same analogy and define a magnonic 
version of this confinement factor based on the probabil¬ 
ity distribution of the magnon confined modes, which, in 
our case is the population inside the guide. 

/■“b 

CF= / IV'nP (6) 

J-Wg 
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FIG. 5. Effect of guide bending on the quality of transport. Subfigure (a)-(c) shows overlaid snapshots of equally spaced 
time-evolved states, where RC is the radius of curvature of the bend and the color axis is (a) Adiabatic bending, in which 
the magnon stays in the transversal ground state throughout the propagation, (b) As we make the bend tighter, the magnon 
starts to excite higher modes, although, it still stays confined in the guide, (c) At small bend radius the magnon leaks out 
of the guide, (d) Population inside waveguide as a function of j/-position, for several radii of curvatures, (e) Bend loss as a 
function of potential depth, emin for several radii of curvature. For a sharper bend, a deeper well is required to avoid bend loss. 


where V'n is the transversal wavefunction of the n-th eigen 
mode of magnon in the presence of the potential. CF 
varies between 0 an 1 and gives a measure of spacial con¬ 
finement. Figure l^a) shows the confinement factor as a 
function of emin for several values of d. A guide with a 
smaller width requires a stronger potential as compared 
to a guide with large width, to achieve the same value of 
CF. Figurej^b) shows the number of confined modes for 
a potential as a function of for several widths. 

As with multi-mode optical fibers, the number of bound 
modes increases with both width and depth of the po¬ 
tential. Along with the spacial landscape, the energy 
landscape also changes with the changing potential. As 
the potential is applied, eigen modes lower their ener¬ 
gies to become bound modes. Figure |^c) shows the first 
ten transversal eigen values as a function of emin- As 
expected, energies decrease with the increasing depth of 
the potential. Figure|^d) shows the energy difference be¬ 
tween the ground state and the first excited state. The 
energy levels also grow apart as they transit from un¬ 
bound to bound modes. 

Similar to an optical waveguide, we can “bend” these 
guides up t o a c ertain bend radius without losing any 
confinemenlP^l^. We show some examples of bending in 
fig-i where we identify three different bending regimes. 
When the radius of curvature of the bend is large, the 


magnon can follow the guide adiabatically, fig. [^a). As 
the radius of curvature reduces, the magnon will leak to 
higher confined modes, fig. [^b). In this kind of bend¬ 
ing, the quantum phase is not preserved, so it is not a 
useful regime for quantum transport. Finally, when the 
bend becomes very sharp (small radius of curvature), the 
magnon cannot follow the guide anymore and it leaks out 
into unconfined modes, fig. [^c). Figure [^d) shows the 
population inside the guide, for several radii of curvature, 
as a function of distance along the propagation direction. 
The bending starts at the 100th site. Radii of curvature 
of 50a and 500a result in partial leakage of the popula¬ 
tion. The radius of curvature of 1000a is a non-adiabatic 
bending in which population oscillates inside the guide 
but does not leak out. Finally the radii of curvatures of 
2000a and 5000a show adiabatic transport. 


IV. MAGNON SPLITTER, AND THE 
MICHELSON INTERFEROMETER 

We now address the issue of the design of a magnonic 
Michelson interferometer. The Michelson interferometer 
comprises an input channel, a splitter element (beam¬ 
splitter or directional coupler), a variable (ideally tune¬ 
able) phase element, and a reflective element that directs 
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FIG. 6. Contours of constant coupling between two spin-guides acting as a directional coupler, as a function of the potential 
depth and the separation between the guides. The coupling strength decreases with increasing depth as well as with increasing 
separation between guides. At small currents, the modes become unbound. The thumbnails on the top and right sides show 
cross sections of the magnetic field perpendicular to the guiding direction, to illustrate the relative change in the shape of the 
potential along both axes. 


the signal back to the splitter. The final output path 
of the magnon then becomes a sensitive function of the 
variable phase. We discuss these elements in turn. 

The directional coupler is an important and commonly 
used device in optical fiber technology. Typically a di¬ 
rectional coupler requires two proximal guides so that 
particles evanescently hop from one guide to the other. 
Similarly, we can achieve the same functionality with our 
magnonic analogs. When two guides are close to each 
other, the magnon can coherently tunnel between the 
guides and this tunnelling rate varies as a function of 
both the distance between the guides and the depth of 
the guides. If the guides are made lithographically, then 
the separation between the guides cannot be varied post¬ 
fabrication, however in our case the depth of the guides 
can be easily modified by varying the magnitude of the 
current in the surface wires. Hence we can realise recon- 
figurable directional couplers. 

Figure shows contours of constant coupling energy 
as a function of current and guide separation. The cou¬ 
pling is stronger when the guides are closer and shallower 
(lower current). To calculate the coupling energy, we per¬ 
formed a full Hamiltonian diagonalisation of a transversal 
cross section of the spin sheet to calculate the transver¬ 
sal ground state of each potential. Then the exchange 
energy was calculated as {L\H\R), where \L) and \R) are 
the ground states of the left and the right guide and H 
is the Hamiltonian with both guides in effect. 


For a Michelson interferometer, we envisage adiabat- 
ically reducing the distance between the guides to ef¬ 
fect the directional coupler, and then adiabatically in¬ 
creasing the distance after half of the population has 
been transferred, i.e. we perform the transformation 
\L) — >• (l/-\/2)(|L) -I- |i?)). The distance that is required 
to transfer half of the population is ^i /2 = Ug7ra/2J^. 

In our model we started the guides at 35 sites apart and 
then adiabatically reduced the distance between them to 
23 sites. The width of each guide was 20 sites, radius of 
the curvature of bend was 5000 sites and Emin was IJ. 
The coupling strength for this particular geometry was 
0.0048J, which gives I 1/2 = 650 sites, fig. Ba). 

Figure shows the evolution where no extra phase 
{(j) = 0) is added to the magnon. Figure [^a) shows 
the propagation up to the reflection, whereas fig. Bb) 
shows the propagation after the reflection. The double 
passing of the beamsplitter has resulted in the magnon 
being shifted from \L) to |i?). 

There are many ways in which a controllable phase 
shift can be introduced. Here we modelled a dynamic 
magnonic crystal-like structures! gf one arm, 

as shown in fig. which provides a tuneable phase shift 
in the right arm. A dynamic magnonic crystal is an ar¬ 
ray of equidistant conducting wires placed on the top of 
magnetic medium. The flow of current in these wires 
causes a periodic magnetic perturbation in the magnetic 
film as shown in fig. (a), which is analogous to a one 
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FIG. 8. (a) The 2 -component of the magnetic field from a 
DMC. DMC was modelled with 10 periodic repetitions with 
inter-wire separation of 20 spin sites. It takes two wires to 
produce one period of the magnetic field, (b) Phase shift as a 
function of current in DMC. When the field is greater 

then J, it acts as a reflecting potential and is no longer useful 
for phase control. Note the break in the horizontal scale. 


FIG. 7. Instances of full time evolution of a Michelson in¬ 
terferometer (for the device schematics see fig. Bb))- White 
lines represent the centre of each guide. Snap shots of the 
population are overlaid on top of each other, where the colour 
represents the I'i/’P- Each snapshot is scaled such that its peak 
value appears as bright as the first instance. The right arm 
has a DMC at its end, which is used as a tuneable source 
of phase difference between both arms (see fig. for de¬ 
tail). Numbers placed next to each snap shot are the time 
stamps (in units of 1/J) and scaling factor respectively, (a) 
The magnon was initialised in state \L) and it splits into 
{1/'/2){\L) -f \R)) upon passing through the splitter. In this 
particular case there is no relative phase shift upon reflection, 
(b) After reflection and second pass through the splitter the 
magnon is transmitted to |i?). 


dimensional magnetic crystal. The depth of the pertur¬ 
bation is controlled by the current in the wire^^. In our 
case the periodicity of the magnonic crystal was chosen 
to be 20 spin sites with 10 periods. As we change the 
current in the crystal, we change the depth of the poten¬ 
tial wells caused by the DMC. The phase change picked 
up by the magnon will be given by the integral over the 
entire DMC potential, which is linear in the minimum of 
the potential, [fig.|^b)]. 

Figure [^b) also shows the reflection and transmission 
of the magnon from the DMC as a function of The 

reflection from the DMC potential will be a function of 
magnon momentum. For a magnon moving with a speed 
2aJ, which is in the desired linear dispersion regime, we 
find for J, there is no appreciable reflection of 

the magnon. As becomes larger than J, it starts 

to reflect the magnon. Note that this reflection threshold 
depends on the speed of the magnon. Slower magnons 


o|i:) + i|i?) i|T) + o|i?.), 



FIG. 9. Population in each arm at the end of the protocol, as 
a function of the depth of the potential due to the magnonic 
crystal, which in turn is a function of current in the crystal. 
The population in right and left arm is shown in red and blue. 
We initialise the magnon in state \L). With no phase shift 
the final state is |J?). As the potential depth increases, the 
population varies sinusoidally between zero and one in both 
arms. The population goes through the whole state space and 
comes back to |L) at /J « 0.07. 


will be reflected from weaker potentials. All of our cal¬ 
culations were done in the limit <C J. 

Reflection was achieved by the hard wall boundary 
condition at the edge of the spin sheet, which reverses 
the direction of the magnon, does not introduce any rel- 
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ative phase between the two arms of the interferometer. 
After reflection, the population in the right arm passes 
the dynamic magnon crystal a second time, doubling the 
effect of its phase shift. For a given device, the final 
state of the interferometer is only a function of the rela¬ 
tive phase between two arms, which in turn is a function 
of current in the crystal. 

Considering the case where the magnon always ini¬ 
tialised in \L), increasing the current in the DMC leads 
to a relative phase shift in the right arm. This leads to a 
sinusoidal variation in the final population in each arm as 
a function of the phase shift, as expected (fig. |^. When 
« 0.017 the final population is equal in each 
arm (l/-\/2)(|L) -|- |i?)). As we keep increasing the cur¬ 
rent, the trend continues and the final state goes through 
|L), (l/-\/2)(|L) — |i?)) and eventually coming back to |i?) 
at / J ~ 0.07. This shows that the current through 
the dynamic magnonic crystal can be used as a tuneable 
element for magnon switching. 


V. REALISTIC SYSTEMS 

Up until now we have presented our results for a 
generic Hamiltonian. To translate our scheme to any 
particular implementation only requires knowledge of the 
coupling strength J, spin separation a and the desired 
wire to spin-sheet separation d. Here we consider a two- 
dimensional spin sheet of nano-patterned phosphorous in 
a silicon latticd^. This system has been widely studied 
due to its relevance to phosphorous in silicon quantum 
computatiori^iJl^ and quantum transporlP^M! 

We take the inter-donor spacing to be a=10nm, for 
which the J-coupling between donors is 40/re\®. If we 
chose the spin guide dimensions Wg = d = 100a = l/rm, 
then according to fig. |^we require emin = 1 x 10“'* J to 
confine a single mode. The current required to produce 
that onsite energy is I/dn = emin, where n = 27r/7fi,/r, 
which gives I = 2.74 ^A. 

For the Michelson interferometer we need two further 
geometrical parameters; the half-coupling length - the 
length that is required to transfer half of the magnon 
population to the other guide - and the adiabatic ra¬ 
dius of curvature, which is the radius at which magnon 
turns without leaking or coupling to the other modes. 
For a high contrast Michelson interferometer we chose 
h /2 = 500a which correspond to an inter-guide separa¬ 
tion of 50a = 0.5^m {I 1/2 = Ug7ra/2J^). For a bending 
guide stronger potential is required to keep the magnon 
confined. If the chosen radius of curvature for interfer¬ 
ometer is 5000a, then the required emin will be 10“^ J 
fig-in: e) which corresponds to a current of 27.4/xA. These 
parameters are well within the experimentally achievable 
limits and would lead to a total device length of 5/iTO. 
There is considerable flexibility in the choice of parame¬ 
ters to build a working device. 

One important parameter to consider is the coherence 
time of the magnon. The magnon coherence time should 


be long enough to allow the magnon to complete the 
return trip through the interferometer with sufficient fi¬ 
delity to observe the interference effects. In the Heisen¬ 
berg framework the maximum speed of the magnon is 
2Ja/H, which gives 607m/s for this particular system. 
This results in 16.7ns round trip time through a device 
of length 5p,m. The most recently reported T-2 in P;Si 
is 2s at 5K temperatur^^, which is orders of magnitude 
longer than the magnon travel time. 

VI. CONCLUSION 

We have proposed a scheme for guiding magnons in 
a two-dimensional ferromagnetic sheet using surface cur¬ 
rent carrying wires. Through numerical simulation based 
on the Heisenberg model, we demonstrated that transver¬ 
sal confinement can be achieved in this setup. We also 
presented a model of a magnonic Michelson interfer¬ 
ometer. The magnon is split and recombined using a 
magnonic equivalent of a directional coupler. The ex¬ 
tra phase was added onto one arm through a dynamical 
magnonic crystal. By changing the amount of current 
in the crystal one can obtain any desired combination of 
population in each arm. Equally, this device could be 
used as a sensor of any field capable of perturbing the 
acquired magnon phase. 

One possibility afforded by our scheme is the design of 
magnonic devices capable of performing non-determinstic 
linear magnonic quantum computation, by analogy with 
non-deterministic linear optical quantum computatioiP^. 
Non-determinstic linear schemes utilise interferometric 
elements and the ‘hidden’ nonlinearity introduced by 
measurement. When magnons are distributed over many 
spins, as we have considered here, than they behave as 
simple type H bosons, and therefore show bunching in 
Hong-Ou-Mandel type configuration^. Therefore our 
results imply that full non-deterministic quantum gate 
operations can be simply ported to the magnonic case. 
One further advantage of our reconfigurable scheme is 
that it is possible to dynamically switch guides on and 
off, and this may lead a natural realisation of schemes 
with ‘shortcuts’ through high dimensional space^^, and 
more generally to non-trivial consideration of optimisa¬ 
tion of the Hilbert-space dimensionality of the resulting 
magnonic circuit for optimal computatiorP^. 
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VII. APPENDICES 


A. Effect of diagonal coupling 



(b) 




FIG. 10. (a) A Square lattice, where the black dots repre¬ 
sent the lattice sites. The black lines is the standard nearest- 
neighbour coupling Je i.e. the coupling between {i,j) and 
(i± l,j) , and between {i,j) and {i,j ± 1). The red lines rep¬ 
resent the next-nearest-neighbour or the diagonal couplings 
Jd, which is the coupling between {i,j) , (i ± 1, j ± 1). (b) 
Dispersion relation of a two-dimensional sheet along the high 
symmetry lines, with (blue-dotted) and without (red) the di¬ 
agonal couplings (Jd). Along the T-X-direction, which corre¬ 
sponds to the principle axis in real space, both lines coincide. 
However, they differ completely along X-M arm and partially 
along the M-F arm of the dispersion relation, (c) A two di¬ 
mensional k-space of a spin-sheet. Two red dots and the arrow 
indicated the initial and final position of a magnon that was 
launched along the principle axis and then bent to 45 degrees. 


Here we will compare the two Hamiltonians by calcu¬ 
lating their dispersion relation through their action on 
an arbitrary wavefunction. A two dimensional Heisen¬ 
berg Hamiltonian with the diagonal coupling is, 


i-1 + f + 1 ~ 

+ V'i+l J-l-1 ~ 4'0ij) 

= (8) 

'ipm,n is the (to, n)th element of the eigne-vector, where 
TO and n corresponds to the spin site indices in x and 
y-direction. 

Using a trial wavefunction '4’{Tn,n) = above 

equation becomes 


^^i[k^{i-a)+kyj] ^i[k^{i+a)+kyj] ^i[k;^i+ky{j-a)] 

_|_ ^ilk^i+ky{j+a)] _ _|_ ;l£^^i[k^{i-a) + ky{j-a)] 

Je 

_|_ ^i[k-i.{i-a)+ky(j+a)] _|_ ^i[k^(i+a)+ky(j-a)\ 

_ P 

_j_ ^i[k^{i+a)Jtky{j+a)] _ 4 ^ _ _ _^i{k:ni+kyj) 

Je. 


_j_ g ikyCL _|_ ^ikyCl _ 4^ 


_j_ ^l^^^—iikxa-\-kya) ^—i{kxa—kya) 

Je 

I ^i^kj.(l kyO,') I Ala; Cl) _ _ 


-E 


( 10 ) 


-E 

—— = 2[cos(fca;a) -I- cos(fcya) -I- 2] 

-\- 2^[cos(A:a;a -I- kyo) + cos{kxa — kyo) -|- 2] (11) 

^ e 


E = 4(Je + Jd) — 2Je[cos{kxa) + cos(A:ya)] 

— 4Jd[cos(fca;o)cos(fcya)] (12) 

If the magnon is only traveling along one direction, say in 
x-direction and ky = 0, then the above equation becomes 


H = -JeE (7) 

'id 

Jd ^ ^ -f .Sj+i j — ij -|- ^ 

hj 

where Je is the standard edge coupling and Jd is the 
diagonal coupling. 

If Ip is the eigen-state of the eq. Q i.e. iJ|'0 >= E\tp >, 
where E is the corresponding eigen-energy, then, expand¬ 
ing the left side of this equation into the action of the 
eq. (|^ on the state ip , 


E =2{Je + 2Jd) [1 - cos(A:,a)] (13) 

This is exactly a dispersion relation of a spin-sheet 
except Je is now Je + 2 J^. 

Figure 10 (a) shows a square lattice with lines rep¬ 
resenting the couplings. Black lines represent the edge 
coupling and red lines represent the diagonal couplings. 
Figure 10 (b) shows the dispersion relation of the spin 
sheet, with (blue-dotted) and without (red) the diagonal 
coupling. Both lines coincide for the F-X part of the 
curve, however the y-axis scales are different for both 
curves. This shows that the addition of Jd acts as 
a scaling factor in the dynamics of a magnon that is 
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traveling along the principle direction in a sheet with 
nearest-neighbour coupling. 


The region of interest in the band structure is the re¬ 
gion that is around the 2J line. This is due to the fact 
that we always launch the magnon with the energy of 2J, 
which corresponds to the maximum gradient and hence 
maximum group velocity, which the magnons travel at. 
Once launched along the principle axis with the wave vec¬ 
tor of 7r/2, which corresponds to the energy of 2J, any 
bending starts to transfer the population along the M — T 
arm of the dispersion curve at fc = tt/S, while keeping 
the energy constant [Fig 10(b,c)]. The M — F arm of the 


dispersion relation, with and without the diagonal cou¬ 
plings, only differs by approximately 3% at the energy of 
2J. Based on this we can approximate the dynamics of 
the magnon in the spin-sheet with a diagonal coupling 
by simply simulating a spin-sheet with a rescaled edge 
coupling. 
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